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The coherent pion production induced by neutrinos in nuclei is studied using a delta hole model 
in local density approximation taking into account the renormalization of A properties in a nuclear 
medium. The pion absorption effects have been included in an eikonal approximation. These effects 
give a large reduction in the total cross section. The numerical results for the total cross section are 
found to be consistent with recent experimental results from K2K and MiniBooNE collaborations 
and other older experiments in the intermediate energy region. 
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The present experiments at K2K and MiniBooNE are 
designed to search for ^-oscillations in disappearance 
and v e appearance channels. In these experiments, the 
i/^-spectrum is determined by the observed energy spec- 
trum of muons which are predominantly produced in for- 
ward direction through the charged current quasielastic 
reactions induced by v^. In this kinematic region the 
major background to these events come from the non 
quasielastic events in which pions are produced through 
coherent and incoherent processes induced by charged 
and neutral weak currents in ^-nucleus interactions. The 
neutral current induced 7r° production is of particular 
importance as it constitutes a major background to the 
electron signal in v e appearance channel. The analysis of 
^-oscillation experiments, therefore, requires a better un- 
derstanding of pion production processes in ^-nucleus in- 
teraction and many experiments are being done to study 
^-induced coherent and incoherent production of pions 
from nuclei 0. 

The coherent production of charged pions from 12 C has 
been recently studied by the K2K collaboration 0. The 
coherent production of neutral pions has been studied 
by the K2K collaboration for 16 O and by the Mini- 
BooNE collaboration for 12 C 0. Many experimental 
groups have earlier studied the coherent production of 
pions induced by neutrinos from other nuclei at higher 
energies 0-0- Theoretically various authors 0-0 have 
used Adler's PCAC theorem to predict the total cross 
sections for neutrino reactions, which overestimate the 
experimental cross sections at low energies. In these cal- 
culations the nuclear medium effects arc included only 
through the final state interaction of the outgoing pions 
with the nucleus using a model of pion nucleus scattering. 
A theoretical framework for treating the nuclear medium 
effects in the neutrino production of coherent pions us- 
ing some model of nuclear structure has been recently 



discussed by some groups |l(| but no definite predictions 
are made for the kinematics of the neutrino oscillation 
experiments of present interest. 

We present in this letter a calculation of the neutrino 
induced production of coherent pions from nuclei at inter- 
mediate neutrino energies. The nuclear medium effects 
are taken into account in the weak production process as 
well as in the final state interaction of the outgoing pions 
with the nucleus. The calculation uses the local density 
approximation to the delta hole model which has been 
used earlier to study photo and electro production of pi- 
ons from nuclei flU ]. The final state interaction of pions 
has been treated in eikonal approximation with the pion 
optical potential described in terms of the self energy of 
a pion in a nuclear medium calculated in this model [l2| . 

In the kinematical region between threshold and few 
hundred MeV excitation energy relevant for the produc- 
tion of pions in neutrino oscillation experiments at inter- 
mediate neutrino energies, the dominant contributions 
are due to the excitation of the A resonance. We treat 
the A as a Rarita Schwinger field '5 \ and use the on shell 
form ^| of the A propagator given by 



Aa ct (P) = 



p8/a (p) 



(1) 



where T is the energy dependent width of A and V^(P) 
is the spin-| projection operator. 
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=0.36 to describe the on 
and the standard model 



shell ANvr interaction [f^-^j 
of weak interactions to describe the A-N isovector charge 
and neutral transition currents corresponding to W and 
Z exchanges. The matrix element for the charge current 
7r + production from the proton, is calculated using the 
Feynman diagrams shown in Fig. 1, with the A propa- 
gator given in eq. (1) and the p — > A ++ weak transition 
current J-^(s) in the s-channel, given by 
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and a similar isospin rotated p — » A transition current 
Jyj{u) for the u-channel diagram. u Q (P) and u(p) are the 
Rarita Schwinger and Dirac spinors for A and nucleon. 
The form factors (j=l-3) and C^(j=l-4) are taken 
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FIG. 1: Feynman diagrams 
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from the experimental analysis of neutrino production of 
pions from proton and deuteron targets The 
q 2 dependence of these form factors has been discussed 
in the literature by many authors 0-0]- We use a 
dipole parametrization for the vector form factors and 
a modified dipole parameterization for the axial vector 
form factors given in Refs. |15| . |17| with a vector dipole 
mass My = 0.73 GeV and an axial vector dipole mass 
M A = 1.05 GeV. 

The transition current for the charged pion production 
from neutron targets is calculated in a similar manner us- 
ing the n — > A + transition current for the s-channel and 
n — > A~ for the u-channel diagrams at the weak vertex. 
For the neutral current induced reaction, we obtain the 
transition current in a similar manner using the p — > A + 
transition current for the production of ir° from protons 
and the n — > A transition current for the production of 
7T° from neutrons in s and u channels. With these con- 



the nuclear transition current J^r^ for the 



siderations 

coherent charged (neutral) pion production from a nu 
cleus A X is obtained by summing over all the occupied 
nucleons in the amplitude and is written as 



J w,z 



M 



•J Pw,zMe s(M -W 



(3) 

Pi is the momentum of the A resonance in i(=s,u) chan- 
nel. T^j z (i) is the nonpole part of the kinematic factors 

involving transition form factors C^' A (g 2 ), pl^ z (r) is the 
linear combination of proton and neutron densities incor- 
porating the isospin factors for charged and neutral pion 
production from proton and neutron targets correspond- 
ing to W and Z exchange diagrams in i(=s,u) channel. 

The mass and width of the A appearing in eq. (3) 
may be modified in the nuclear medium due to inhibi- 
tion of the A decay caused by the Pauli blocking of nu- 
cleons in the A — >N7r decay and due to two and three 
body A absorption processes like AN— >NN, AN— > AN-7T 
and ANN— >NNN, which may take place in the nuclear 
medium. These modifications in the A properties have 
been studied by many authors 01 1 E3 by calculating 
the self energy £a of A in the nuclear medium leading 
to density dependent changes in Ma and T, given by 



FIG. 2: cr(E„) vs E„ for coherent tt + production in 12 C (see 
text for details). 
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FIG. 3: o-(E„) vs E„ for coherent n° production in 12 C(solid), 
27 Al(dashed) and Freon(dotted) with nuclear medium and 
pion absorption effects, along with the experimental re- 
sults 0-13 



M A = M A + Re£ A and ^ = ^ - ImS A) where T is 
the A width modified due to Pauli blocking. We use the 
expressions of Oset and Salcedo 0] for T, ReSA and 
ItoEa where the density dependence of Ea and T in the 
nuclear medium has been studied in a local density ap- 
proximation to the delta hole model. Thus incorporating 
the nuclear medium effects, the A-dependent hadronic 
factors in eq. (3) become density dependent and the nu- 
clear transition current z is now modified to 
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The final state interaction of the pion is taken into 
account by replacing the plane wave pion by a distorted 
wave pion. The distortion of the pion is calculated in the 
eikonal approximation in which the distorted pion wave 
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FIG. 4: jZj vs - q 2 at E„ = l GeV, for coherent n + pro- 
duction in 12 C nucleus without (dotted), with(dashed) nuclear 
medium effects and with nuclear medium and pion absorption 
effects(solid). In the inset the final results for < -^j vs — q 2 
averaged over the K2K and MiniBooNE spectrum are shown. 
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FIG. 5: -jZ- vs pTv at E v =l GeV. The various curves corre- 
spond to the cases described in Fig. 4 



function is written as 



'■■5^ 25bn( P (M'))<fe' 



(5) 



where f = (6, z). II(p) is the self energy of the pion 
calculated in the local density approximation of the delta 
hole model and is taken from Ref. hj 



U(p) = 



4 f Una \ 2 M 2 



W- M A + 



(6) 



Note that we use the nonrelativistic form of the energy 
denominator in the delta propagator for calculating the 
pion self energy. 

The nuclear transition current operators z dis- 
cussed above in eq. (3), with incorporation of nuclear 
medium effects through eq. (4) and the final state interac- 



tion effects through eq. (5), are used to numerically eval- 
uate the total cross sections er(E„), differential cross sec- 
for 12 C, 27 Al and Freon(CF 3 Br). We 



tions 



take the proton density p p (r) = \p(r) and the neutron 
density p-air) = ^j^-p{r) which should be appropriate for 
describing nuclcon densities in intermediate mass nuclei. 
The nuclear densities p{r) are taken from Ref. [20| . 

In Fig. 2, we have presented the results for the total 
scattering cross section (cr cc ) for the coherent charged 
current reaction induced by in 12 C. The results for 
cr (E„) (scaled by a factor of i) vs E y are shown with- 
out nuclear medium effects and with nuclear medium ef- 
fects. When the pion absorption and nuclear medium ef- 
fects, are both taken into account the results for a cc (E v ) 
are shown by the solid lines. We see that the nuclear 
medium effects lead to a reduction of 30-35% around 
E„=l-1. 5 GeV in (T CC (E„) while the reduction due to 
final state interaction is quite large. We have also shown 
in this figure the results for cr cc (E„) vs E„ when a 
cut of 450 MeV is applied on the muon momentum as 
done in the K2K experiment. We obtain a total cross 
section of 6.93x 10 _40 cm 2 in 12 C which corresponds to 



0.578x10 



-40 



nucleon • 



and is consistent with the exper- 



imental result of cr cc <0.642x 10" 40 - 
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ported by the K2K collaboration 2]. Our calculations 
show that: 

(i) the contribution to the cross section comes mainly 
from the s-channel diagram(> 90%) which is dominated 
by on shell A, thus making the off shell correction quite 
small. The inclusion of off shell effects by introducing a 
form factor at the 7T.ZVA vertex 0-01 leads to a reduc- 
tion in the cross section which is estimated to be 4-6 % 
in the energy region of 1-2 GeV. 

(ii) the contribution to the cross section from the vector 
current is negligibly small (<2 %) and the major contri- 
bution comes from the axial current only, leading to near 
equality of neutrino and antineutrino cross sections. 

We show in Fig. 3, the total cross section ct nc (E j/ ) 
for the neutral current induced ir° production from 12 C, 
27 A1 and CF3Br(Freon), along with the experimental 
results from the MiniBooNE collaboration for 12 C jj], 
from the Aachen collaboration for 27 A1 [j| and from the 
Gargamelle collaboration for Freon We see that the 
theoretical results for the neutral current induced co- 
herent 7T° production are in reasonable agreement with 
presently available experimental results in the intermedi- 
ate energy region. 

We also present in Figs. 4 and 5 the differential cross 
sections ^ and for charged pion production at 
E„=l GeV where nuclear medium and final state interac- 
tions effects arc shown explicitly. In the inset we exhibit 
the differential cross sections < ^ > and < j^- > av- 
eraged over the K2K and MiniBooNE neutrino spectra 
without applying any cuts on the muon momentum. 

We would like to emphasize the important role of nu- 
clear medium effects which reduce the total cross sections 
and help to obtain reasonable agreement with the ex- 
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perimental values. The uncertainty in calculating these 
nuclear effects conies mainly from uncertainties in the 
value of the A-self energy which is calculated [l2j within 
an accuracy of 15-20 %. This leads to an uncertainty of 
8 — 12% in the total cross section. There is an additional 
uncertainty due to various parameterizations used in the 
literature for the weak transition form factors |18| , which 
is small (3-5 %) as this process is dominated by differen- 
tial cross sections at very low <? 2 (Fig. 4). Thus the total 
uncertainty in the present calculation is estimated to be 
about 15%. 

Similar results for a ' NC (E v ) in 12 C have been re- 
cently obtained by Paschos and Kartavtsev [jj who have 
calculated only the axial vector contributions to the total 
cross sections. They use Adler's PCAC theorem 8] to re- 
late the forward neutrino cross section to the 7r- 12 C elas- 
tic scattering and extend it to non-zero Q 2 (< 0.2 GeV 2 ). 
The experimental data on ir- 12 C elastic scattering cross 
sections have been used, which are available only in a 
limited region of pion energy and momentum transfer. 
Moreover, for other nuclei where no experimental cross 
sections are available an approximate scaling law for the 
A-dependence of pion-nucleus elastic scattering cross sec- 
tions has been used to obtain cross sections for these 
nuclei from 12 C data. In their model no nuclear ef- 
fects are included in the weak pion production process 
and the final state interaction of pions is included only 



through the dynamics of pion- 12 C elastic scattering. On 
the other hand, in the present work the results obtained 
from a microscopic study of the A-nucleus and pion- 
nucleus dynamics have been used to evaluate various nu- 
clear medium effects in the weak pion production process 
as well as in the final state interactions of pions with a 
nucleus. 

To summarize, we have studied in this letter total cross 
section a(E u ), differential cross sections ^ and for 
the neutrino induced coherent production of charged and 
neutral pions in a model, which takes into account nu- 
clear medium effects in the weak pion production process 
through A dominance treated in local density approxima- 
tion. The final state interaction of pions with the nucleus 
is described in an eikonal approximation with a pion op- 
tical potential derived in terms of the pion self energy in 
the nuclear medium. The good agreement between the 
theoretical and experimental results in the energy region 
of 1 GeV is obtained mainly due to inclusion of nuclear 
medium and pion absorption effects. The method may be 
useful to analyze the neutrino induced pion production 
data at neutrino energies relevant for neutrino oscilla- 
tion experiments being done by K2K, MiniBooNE and 
J-PARC collaborations. 
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